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ABSTRACT 

Cosmological simulations of magnetic fields in galaxy clusters show that remarkable agreement between simula- 
tions and observations of Faraday rotation and radio haloes can be achieved assuming that seed fields of ~ 10~ 9 G 
were present at redshifts ~ 15-20. The structure of the seed field is irrelevant for the final intracluster field. On 
average, the field grows exponentially with decreasing cluster redshift, but merger events cause steep transient in- 
creases in the field strength. Typical field-reversal scales are of order 50/z kpc. In most cases, the intracluster 
fields are dynamically unimportant. Assuming secondary electron models, the average structure of cluster radio 
haloes can naturally be reproduced. 



1. Introduction 



Magnetised plasmas are optically birefringent. The orientation 
r—i of polarised light traversing a magnetised plasma is therefore 
■^j" changed; this effect is called Faraday rotation. Being propor- 
tional to the square of the wavelength, it can be detected by ob- 
servations of polarised radio sources in two or more frequency 
CD bands. Numerous such observations show that clusters contain 
magnetic fields which are smooth, i.e. ordered on scales com- 
parable to the cluster scale (e.g. Vallee, MacLeod & Broten 
^ 1986, 1987; Vallee 1990; Dreher et al. 1987; Kronberg 1987; 
Q^im et al. 1990; Kim, Kronberg & Tribble 1991; Clarke, Kro- 
i nberg & Bohringer 2001 to name only a few). Micro-Gauss 
O field strengths are routinely inferred. Figure |l] illustrates data 
-4— » from Clarke et al. (2001). The Faraday rotation measure seen 
£2 in polarised radio sources is plotted against their projected sep- 
. . aration from the nearest Abell cluster. Evidently, the rms ro- 
J> tation measure increases substantially towards cluster centres, 
k> wni le it is compatible with zero at cluster-centric distances be- 
rN yond ~ 1 /i _1 Mpc. 

This demonstrates that at least part of the Faraday rotation 
measured in sources behind galaxy clusters is due to intracluster 
magnetic fields rather than magnetic fields which are intrinsic to 
the sources themselves. 

The Faraday rotation measure is proportional to the integral 
along the line-of-sight over the parallel component of the mag- 
netic field, times the electron density. The same rotation mea- 
sure can therefore be caused by strong fields which often change 
their direction, thus causing substantial cancellation, or weak 
fields which are almost homogeneous across a cluster. It is thus 
far from unique to infer magnetic field strengths from Faraday- 
rotation measures. Typically, fields with strengths of order /jG 
are inferred, assuming simple magnetic -field geometries and 
typical electron densities derived from cluster X-ray emission. 

Another manifestation of intracluster magnetic fields is pro- 
vided by radio- and hard X-ray emission. Radio emission is 
caused by relativistic electrons gyrating in the intracluster field. 
The same population of relativistic electrons must Compton- 
upscatter photons of the cosmic microwave background into the 
hard X-ray regime. From the absence, or the low intensity, of 
this non-thermal, hard X-ray emission, lower limits can be in- 
ferred on the intracluster magnetic fields, because if the rela- 
tivistic electron population is poor, the magnetic field has to 
be strong in order to explain the cluster radio haloes. These 
lower limits are typically also of micro-Gauss order (e.g. Baz- 
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Fig. 1. — Data by Clarke et al. (2001) showing the Faraday rotation 
measure in radio sources against their projected distance from the near- 
est Abell cluster. While the rms rotation measure is compatible with 
zero beyond ~ 1 Mpc, it increases substantially towards smaller cluster- 
centric distances. This demonstrates the existence of smooth intraclus- 
ter magnetic fields. 



zano 1990; Fusco-Femiano et al. 1999; Bagchi et al. 1998; Hen- 
riksen 1999). 

This outlines the motivation for our numerical work sum- 
marised below. We wish to understand how magnetic fields 
evolve in galaxy clusters, how they act back on their host clus- 
ters themselves, how the initial conditions need to be set up so 
that observables like rotation measures and radio haloes can be 
reproduced, and how the magnetic fields are structured inside 
clusters. 



2. Numerical Technique 

For these simulations, we use the GrapeMSPH code (Dolag, 
Bartelmann & Lesch 1999), developed from M. Steinmetz' 
(1996) GrapeSPH code. The code uses the Grape hardware 
(Ito et al. 1993) to trace dark-matter particles according to the 
softened force law 

where a is the index of the particle under consideration, r a and 
v fl are its position and velocity, and e„ is its softening length. 
The summation extends over all other particles at positions r, 
with masses mi. 

The dynamics of the gas is treated in the SPH approxima- 
tion (Lucy 1977; Monaghan 1992). Particles are replaced by 
extended, soft spheres whose density is described by a kernel 
function W with variable width h. Any physical quantity A at 
position r is approximated as 

(A^^m^Wir-r^h), (2) 
i Pvi) 

where the sum extends over all particles which come sufficiently 
close to r. The code combines the advantages of SPH, which 
is a fast method with adaptive resolution, with those of the 
Grape hardware, which provides supercomputer performance 
on low-end hardware. In our case, the host workstation is a two- 
processor SUN Ultra-Sparc 2300 equipped with five Grape-2>M 
boards. 

The equation of motion for the gas particles can symbolically 
be written 

dv a 

= pressure gradient + viscous stresses + magnetic force 

dt 

+ gravity + D. A . (3) 

In particular, the back-reaction of the magnetic field on the gas 
flow through the Lorentz force is included in the calculation. 
The term £2a symbolises the additional acceleration in presence 
of the cosmological constant. In addition to the momentum con- 
servation expressed by Eq. ([}]), the energy equation is solved. 
The gas is treated as an ideal gas with an adiabatic index of 5/3. 
Cooling and heating are foreseen in the code, but have so far 
been used for testing purposes only. 

We assume that the conductivity of the intracluster plasma is 
infinite, so that the induction equation of ideal magnetohydro- 
dynamics can be used, 

dB 

— = Vx(vxB). (4) 
at 

Equation (Q) and the inclusion of the Lorentz force in Eq. (||) 
constitute our main complements to the original GrapeSPH 
code. 

The numerical viscosity is chosen such as to minimise viscous 
angular momentum transport in particle encounters (Monaghan 
& Gingold 1983; Balsara 1995; Steinmetz 1996). The viscosity 
vanishes exactly in pure shear flows. Finally, all SPH expres- 
sions are symmetrised against particle exchange. 

3. Simulations 

We tested the code extensively, as described in detail in Dolag et 
al. (1999). Specifically, the code succeeds in computing spheri- 
cal collapse models, and in solving shock tube problems (Brio & 



Wu 1988). Although the code does not by construction conserve 
V • B, it turns out in our simulations that V • B is always many or- 
ders of magnitude smaller than \B\ divided by the correlation 
length of B. 

We perform simulations in two different cosmological mod- 
els, SCDM with D.q = 1 and £2 A = 0, and ACDM with Q. = 0.3 
and £2a = 0.7. In both cases, we initialise dark matter dis- 
tributions from the COBE-normalised CDM power spectrum. 
Galaxy clusters are surrounded a layer of more massive parti- 
cles in order to mimic the gravitational tidal field exerted by 
neighbouring large-scale structure. The total particle number is 
limited by the Grape architecture, but our simulated clusters are 
generally resolved by > 10 4 particles. 

Initial conditions are set up such that the average magnetic 
field strength at high redshifts (15 to 20 depending on cosmol- 
ogy) are of nano-Gauss order. In lack of any detailed predic- 
tion of the magnetic field structure at such redshifts, we use two 
vastly different assumptions. The field is either assumed to be 
homogeneous across the entire simulation volume, or chaotic, 
in which case its amplitude is drawn from a pre-defined power 
spectrum, and its direction is random, subject only to the condi- 
tion that V -B = 0. Figure H shows one evolutionary stage of one 
of our cluster simulations as an example. 




Fig. 2. — One evolutionary stage of one of our cluster simulations. 
Three boxes are shown which zoom into the cluster region from left to 
right. The right-most box has a side length of 5 h~ 1 Mpc. It shows a grey 
iso-density surface of the intracluster gas distribution in the centre; the 
grey contours on the box sides are isodensity curves of the intracluster 
gas. The coloured patches on the box sides show the Faraday rotation 
measure, with values ranging within ±500radm~ 2 . Red patches have 
negative, blue patches positive rotation measure. The purple frame in 
the middle of the box shows a slice through the intracluster magnetic 
field, with the arrows indicating field strength and orientation. 

The Faraday rotation measure created by the intracluster mag- 
netic field is show as the coloured patches projected onto the 
box sides. The rotation measures range within ±500radm~ 2 , 
increasing from red to blue. Obviously, the rotation measures 
reach values typical for the observations, they cover a substan- 
tial portion of the cluster's cross section, and they show structure 
on scales of order ~ 100/j^'kpc. The figure also shows how the 
cluster is embedded into the surrounding large-scale structure, 
and it displays the strength and orientation of the magnetic field 
in a slice cut through the cluster core. 



4. Results 



4.1. Field Structure and Amplification 

Starting from nano-Gauss fields at redshifts between 15 and 20, 
intracluster fields reach micro-Gauss strength in and near clus- 
ter centres at low redshift. Such field amplifications exceed ex- 
pectations from simple spherical collapse models of a gas cloud 
with frozen-in magnetic field by about an order of magnitude. 
This becomes possible due to shear flows in the intracluster 
gas, which stretch, bend and entangle magnetic field lines and 
thereby increase the magnetic field strength. The left panel in 
Fig. |^ shows that the magnetic field strength, averaged within 
the central regions of many clusters in our simulated sample, 
increases approximately exponentially with decreasing redshift, 
(\B\) oc \Q- 2 < (Dolag et al. 1999; Dolag, Bartelmann & Lesch 
2001). 

We find that the initial field structure is entirely irrelevant for 
the final field structure after cluster collapse. Simulations start- 
ing from the homogeneous magnetic field set-up lead to final 
magnetic fields which are statistically identical to those obtained 
from chaotic initial fields. This shows that the final structure of 
intracluster fields is completely determined by the cluster col- 
lapse, in the course of which any information on the initial field 
structure is erased. 

The autocorrelation function of the magnetic field reveals typ- 
ical field-reversal scales of order 50/i _1 kpc. The right panel 
in Fig. H shows the power spectrum of the modulus of the fi- 
nal magnetic field strength for various initial conditions. The 
power spectrum resembles a power law for values of k suffi- 
ciently larger than the Nyquist frequency of the box, quite inde- 
pendent of the initial field set-up. The power-law is fairly steep, 
with exponents ranging within —2.5 to —3. This is steeper than 
expected for Kolmogorov turbulence, but compatible with two- 
dimensional magnetohydrodynamic turbulence. It cannot be ex- 
cluded at this stage that the numerical viscosity of the SPH code 
contributes to steepening the magnetic-field power spectrum. 

4.2. Faraday Rotation 

Simulations in both cosmologies reproduce the statistics of ob- 
served Faraday-rotation measures very well. On the whole, 
the amplitudes and coherence lengths of the simulated rotation- 
measure maps, and their projected radial profiles, agree excel- 
lently with the available data. As noticed before, the initial field 
configuration is irrelevant in that respect. Setting up fields with 
nano-Gauss strength at redshifts 15-20 is sufficient. 

Merger events in the cluster history are prominently reflected 
in the spatial distribution of the Faraday rotation measure. To 
illustrate this point, we show in Fig. |^ the fraction J of the to- 
tal cluster cross section which is covered by Faraday rotation 
measures with absolute values exceeding given thresholds. The 
total cluster cross section is defined as the area emitting 90% 
of the total cluster X-ray luminosity. In other words, we study 
the measure of the excursion set of the Faraday rotation measure 
map. 

Figure ^ shows a steady increase in J with decreasing red- 
shift, underlying steep, transient increases which are indepen- 
dent of the threshold of the excursion set. Correlating these in- 
creases with the cluster evolution shows that J increases during 
merger events, reaching its maximum during core passage of the 
merged substructure, and decreasing thereafter as the cluster re- 
laxes. Merger events therefore leave a pronounced imprint on 
the appearance of a cluster in Faraday-rotation measurements. 

Observationally, the rms rotation measure observed in galaxy 
clusters is closely correlated with their X-ray surface brightness. 
While this correlation contradicts simple analytic models for the 
intracluster magnetic fields, it is naturally reproduced in our sim- 
ulated clusters (Dolag et al. 2001). 



BigA 

1 .2 [ 1 1 1 i 1 1 1 i 




0.0 0.2 0.4 0.6 0.8 



FIG. 4. — Fraction J of the total cluster cross section covered by Fara- 
day rotation measures whose absolute values exceed the thresholds in- 
dicated in the plot. The total cluster cross section is defined to be the 
area emitting 90% of the total cluster luminosity. Irrespective of the 
threshold, the curves show steep, transient increases of J on top of 
a more gentle increase. The maxima coincide with the core passages 
following merger events. 



4.3. Dynamics 

The intracluster magnetic fields provide pressure support adding 
to the thermal gas pressure. Hydrostatic equilibrium with the 
cluster potential wells can therefore be achieved with lower gas 
temperatures. However, the fields in our simulated clusters are 
typically so weak that the temperature reduction is of order a 
few per cent only. Although it thereby increases the scatter in 
the mass-temperature relation to some degree, this effect is gen- 
erally unimportant (Dolag, Evrard & Bartelmann 2001). 

Almost everywhere in the simulated clusters, and almost al- 
ways during their evolution, the magnetic fields are dynami- 
cally unimportant compared to the uncertainties in the standard 
P model commonly used to infer cluster masses under the as- 
sumption of hydrostatic equilibrium. Lower cluster tempera- 
tures cause p-model mass estimates to be biased low. Figure || 
shows that cluster masses inferred from p models to the X-ray 
emission of magnetised galaxy clusters can be as low as ~ 60% 
of the true cluster masses, but only in cluster cores and if the 
fields are comparatively strong. However, masses can substan- 
tially be underestimated in highly substructured or merging clus- 
ters (Dolag & Schindler 2000). 

4.4. Radio Haloes 

A natural consequence of intracluster magnetic fields is the de- 
velopment of smooth radio emission on cluster scales, so-called 
cluster radio haloes. Relativistic electrons in the intracluster 
plasma gyrate in the magnetic field and emit synchrotron ra- 
diation. Various models exist for the origin of the relativistic 
electron population. They can be primarily accelerated, but it 
is more likely that they from protons which are accelerated to 
relativistic energies in intracluster shocks. These protons decay 




Fig. 3. — Left panel: Growth with redshift of the magnetic field strength averaged within many of our simulated clusters. On average, the field 
strength increases exponentially with decreasing redshift, approximately <x 10~ 2j . - Right panel: Power spectra for the modulus of the intracluster 
magnetic field strength for various initial field set-ups. Above the Nyquist frequency of the simulation box, the power spectra are fairly steep power 
laws with exponents between —2.5 and —3. 



into relativistic pions and electrons. Assuming such a secondary 
electron model, and putting between 5% and 15% of the ther- 
mal energy content into relativistic protons, our simulations can 
straightforwardly reproduce the radio halo of the Coma cluster 
(see Fig. ||). 

The figure shows the radial profile of the radio surface bright- 
ness in one of our simulated clusters (red curves) and in the 
Coma cluster (blue curve). While Coma's radio-halo profile is 
well reproduced by one particular simulation, the mean profile 
averaged across our complete simulated cluster sample is some- 
what steeper. However, this can easily be modified by assuming 
that different spatial distributions of the relativistic protons from 
which the secondary electrons originate. 

Remarkably, though, the very steep correlation between the 
cluster radio power and the X-ray temperature is naturally and 
accurately reproduced by our simulations (Dolag & EnBlin 
2000). 

5. Summary and Discussion 

We have performed cosmological simulations of magnetised 
galaxy clusters in two CDM cosmogonies, SCDM and ACDM. 
The Grape hardware allowed these simulations to be efficiently 
executed on low-end hardware. We developed our GrapeM- 
SPH code (Dolag et al. 1999) by augmenting Steinmetz' (1996) 
GrapeSPH code with the induction equation of ideal magneto- 
hydrodynamics, and the Lorentz force for the back-reaction of 
the field on the intracluster gas. The code successfully passed an 
extensive suite of test runs. 

Our cluster simulations lead to the following results: 

• Micro-Gauss fields at low redshift are reached starting from 
nano-Gauss fields at redshifts 15-20, depending on the 
cosmological model. Such field amplifications exceed by 
about an order of magnitude expectations from magnetic 



flux conservation in the simple spherical collapse model. 
They are possible due to shear flows in the intracluster gas. 

• The initial structure of the magnetic fields is entirely unim- 
portant for the final structure. Simulations starting from 
homogeneous magnetic fields lead to final fields which are 
statistically identical to those obtained from chaotic initial 
fields. The process of cluster formation wipes out any in- 
formation on the initial field configuration. 

• On average over many independent cluster simulations, the 
magnetic field strength grows exponentially with decreas- 
ing redshift. 

• The simulations reproduce the observed statistics of 
Faraday-rotation measures very well. Cluster merger 
events lead to steep and transient increases in the area of 
the cluster which is covered by substantial rotation mea- 
sures. This area is maximised during core passage of the 
merged substructure. 

• Almost everywhere in the simulated clusters, and almost 
always during their evolution, magnetic fields are dynam- 
ically unimportant. They cause some increased scatter in 
the relation between X-ray temperature and mass. During 
merger events, however, cluster masses inferred from p fits 
can be lower than the true masses by factors of two or more. 

• Assuming a secondary model for the production of rel- 
ativistic electrons, and putting 5%-15% of the thermal 
cluster energy into relativistic protons, our simulated clus- 
ters have radio haloes which agree well with observed 
haloes. In particular, the very steep correlation between ra- 
dio power and X-ray temperature is naturally reproduced. 

On the whole, it seems that our magnetohydrodynamic, cos- 
mological cluster simulations capture many features displayed 




FIG. 5. — Masses inferred from (3 models fit to magnetised clusters, di- 
vided by the true cluster masses, as a function of cluster-centric radius. 
For clusters with comparatively strong fields, the inferred mass near the 
cluster centres can drop to ~ 60% of its true value. Generally, however, 
the uncertainty is small compared to the uncertainty in the (3 model. 



FIG. 6. — Radial profiles of the radio surface brightness for the Coma 
cluster (blue curve) and one of our simulated clusters assuming a sec- 
ondary electron model (red curves). While this particular simulation 
reproduces Coma's radio-halo profile very well, the simulated radio 
profiles are on average somewhat steeper. 



by observed galaxy clusters. However, the resolution of the sim- 
ulations needs to be increased before detailed and reliable pre- 
dictions can be made on the small-scale structure of intraclus- 
ter magnetic fields and on the structure and emissivity of radio 
haloes. At a later stage, cooling of the intracluster gas has to 
be included, because it is expected to change substantially the 
structure and dynamical importance of magnetic fields in cluster 
cores. 
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